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ABSTRACT 

Near-IR observations are important for the detection and characterization of exoplanets using the 
transit technique, either in surveys of large numbers of stars or for follow-up spectroscopic observa- 
tions of individual planets. In a controlled laboratory experiment, we imaged ~ 10 4 critically sampled 
spots onto an Teledyne Hawaii-2RG (H2RG) detector to emulate an idealized star-field. We obtained 
time-series photometry of up to ~ 24 hr duration for ensembles of ~ 10 3 pseudo-stars. After rejecting 
correlated temporal noise caused by various disturbances, we measured a photometric performance of 
<50 ppm-hr -1 / 2 limited only by the incident photon rate. After several hours we achieve a photon- 
noise limited precision level of 10 ~ 20 ppm after averaging many independent measurements. We 
conclude that IR detectors such as the H2RG can make the precision measurements needed to de- 
tect the transits of terrestrial planets or detect faint atomic or molecular spectral features in the 
atmospheres of transiting extrasolar planets. 

Subject headings: instrumentation: detectors, near-infrared - methods: laboratory - methods: data- 
analysis - techniques: photometric - planets:transits 



1. INTRODUCTION 
The dramatic results of the Kepler and COROT mis- 



2011b). Spitzer and HST photometry and spectroscopy 



201l| |Baglin et al.||2009[ ) as well 



sions (Borucki & Koch 

as earlier ground- based surveys have demonstrated the 
power of the transit method for discovering and char- 
acterizing other planetary systems. The combination of 
a planetary radius from a transit measurement with a 
mass from a radial velocity measurement yields a planet's 
density and thus clues to its composition and evolution- 
ary history. Optical and infrared observations can de- 
termine a planet's albedo and equilibrium temperature 
and even probe global atmospheric circulation patterns. 
Spectroscopic transit observations have identified atomic 
and molecular species, including H 2 0, CO2, and CH 4 in 
the atmospheres of planets orbiting nearby bright stars. 

The critical instrumental requirement for all these ob- 
servations is highly precise differential photometry since 
the transit of a Jupiter in front of a solar type star (1 R ) 
causes only a 1% decrease in brightness while the tran- 
sit of an Earth in front of the same star causes only a 
0.0084% dip (where 0.01% = 100 ppm). Wide-field sur- 
veys from the ground are typically limited in precision 
by residual at mospheric effects to 2 ,000^-4,000 pp m, e.g. 
Super- WASP flPollacco et aTj2006 ) and HATNet ( |Bakos| 
et al.||2009[ ). Higher precision measurements, < 1,000 
ppm, has been demonstrat ed with concentrated obser- 
vations of individual stars (Irwin et al. 2010| [Johnson 
et al.||20Q9 ). From space, however, much higher levels ol 
precision can be achieved for hundreds of thousands of 
objects simultaneously. The Kepler and COROT mis- 
sions have demonstrated photometric precision in the 
visible using CCD arrays as low as 15 p pm, limited in 
many cases by stellar variability itself (Jenkins et al. 



using earlier generations of near- and mid-IR detectors 
have demo nstrated prec ision around 100 ppm (Knutson 
et al. 120071 ISwain et al.|2009( iDeming et al.||20(Jt>|). 



There are many advantages to carrying out transit ob- 
servations in the near-IR compared to the visible band: 
the photon fluxes from K and M stars are higher in the 
near-IR, up to a factor of 50-100 for the coolest M stars; 
total photon fluxes from hotter F and G stars are fac- 
tor of ~2 higher in a wavelength band that extends out 
to 1.6 /im rather than in a visible band that cuts off at 
the 1 /im CCD limit; starspot noise is reduced compared 
to visible wavelengths due to the lower contrast between 
areas with typical s pot and photospheric temperatures 
(Frasca et al.||2009|); and the effects of limb- darkening 



are greatly reduced, simplifying e xtraction of planeta ry 
parameters from the observations ( |Desert et al.||201l| ). 

The detection of transiting K and M stars is of great 
scientific interest since the smaller radii of late type stars 
en hance s the signals from Earth-sized planets (Gould et 
aL||2003 ). The transit of an Earth in front of an M8 star 



produces a 1% brightness dip, 100 times larger than for 
a G2 star. Additionally, the location of the habitable 
zone (HZ) where a planet's temperature might allow for 
the existence of liquid water, moves in closer to the star, 

anz oc L l J 2 oc M^ 2 with a ~ 3.4, increasing both the 
probability of a geometrical alignment and the frequency 
of transits (a few days or weeks compared to a year for a 
G2 star). Initial results from the Kepler mission suggest 
that low mass stars have an increased incidence of low 



mass planets ([Howard et al.|20TT ) further enhancing the 
importance of later spectral types. Figure [I] highlights 
the increased detectability of transiting planets in the 



2 



100,000 




GO KO K5 M2 M4 M5 M6 

Spectral Type 



Figure 1. The 10-1,000-fold Cool Star Advantage of lower luminosity K and M stars relative to a GO star for detecting habitable zone 
(HZ) transits is illustrated by a geometrical metric which takes into account the increased alignment probability, the increased frequency 
of transits, and the increased transit depth as a function of spectral type (black dashed line; |Gould et al.] ( |2003| >). The Near-IR Advantage 
(blue dotted line) is the ratio of the number of stars detectable at a given signal to noise ratio (SINK) m a broad J\y passband (0.66-1.70 
^m) with H2RG detectors compared to the number detectable with CCDs working in the visible (0.6-1.0 /im). A near-IR survey could 
detect twice as many planets orbiting F and G stars and between 3-7 times more planets around K and M stars. Finally, a Monte Carlo 
simulation discussed in §5 shows that ~115,000 stars in an all-sky survey lasting 1.5 years would be bright enough for the detection of an 
1.5 planet located in the habitable zone (HZ) of its host star (solid red curve). In this context a detection requires the observation of 
at least 3 transits and a mission total SNR>7. 



habitable zones of late type stars. The recognition that 
cool stars offer an excellent opportunity for detecting HZ 
planets has lead to ground-based searches directed at M- 
dwarfs ( |Charbonneau et al.||2009 ) and generated consid- 
erable interest m improving near-IR Doppler precision 



Bean et al. 2010 Plavchan et al. 



Rein ers et al. 2010 

|2oTT|). 

The large signal from a planet transiting a cool star 
also enhances the spectros copic potential of the se sys- 
tems. As pointed out by |Deming et aL (2009), tran- 
sit observations using the James Webb Space Telescope 
(JWST) could identify numerous atomic and molecular 
species in the atmosphere of a ~2 R e planet transiting 
a M5 star. Future missions have also been proposed to 
exploit the potential of spectroscopy to characterize exo- 
planet atmospheres. A critical need for successful transit 
spectroscopy are planets orbiting nearby, bright stars so 
that the necessary precision can be achieved in a rea- 
sonable time. Kepler and COROT targets are too faint 
for this purpose and planets discovered by ground-based 
transit surveys are too few in number and limited to 
relatively large radii due to the constraints of observing 
through the atmosphere. 

Because nearby bright stars are distributed isotropi- 



cally, there is interest in conducting an extremely wide 
field-of-vi ew panoramic survey using a low cost orbital 

Beich- 

man et al 



platform ([Ricker et al.||2010| |Deming et al.||2009 
a survey would target 1-2 



illion 

stars which would be at least 5 magnitudes brighter than 
the Kepler stars. Such a survey could be conducted with 
multiple 12.5 cm apertures, a factor 10 smaller than Ke- 
pler's 1 m telescope, o perating either in the visible us- 
ing CCD focal pla nes ( | Ricker et al. [20 10 Deming et al. 
2009; jClaudi||2010D or in the near mtrT 

M 



al 



Claudi|[2010| or in the near infrared! Beichman et 
ilting the observational band from the vis- 



ible (0.6-1.0 jam) into the near-IR (0.66-1.65 /imj^] for 
otherwise identical survey facilities results in a 3 to 10- 
fold increase in the numbers of cool stars (K5-M5) bright 
enough for the detection of ~ 1 Rf; planets (dashed 
line in Figure [TJ. A Monte Carlo simulation (§5) shows 
that an all-sky survey of 2 million stars would include 
~ 150, 000 stars bright enough to allow the detection of 
planets as small as 1.5 Rq located in the habitable zone 



1 The short wavelength cutoff is selected to exclude Ha emission 
which can be variable in cool stars. The long wavelength cutoff is 
tailored to exclude small temporal or spatial variations in detector 
QE at the limit of its responsivity. 
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(HZ) of its host star with a mission total signal to noise 
ratio of 7 or higher. The vast majority of these stars have 
spectral types later than MO and are thus preferentially 
detected in the near-IR. 

The increased near-IR yield over the visible comes from 
the increased signal-to-noise ratio (SNR) in the near- 
IR due to st ronger photo spheric emission in the near- 
IR ( |Hauschildt et al.|1999 ) and/or broader measurement 



bandwidth. The star counts above a given flux threshold 
increase as the photon luminosity in a bandpass (La A A) 
to the 3/2 power. The number of stars above a given 
SNR threshold goes as TV oc S7ViT 3 (L A AA) 3 / 2 , and thus 
by increasing La A A, the number of stars above a fixed 
SNR increases as (L\A\) 3 / 2 . Figure [l] also shows that 
the near-IR measurements do as well or better than a 
comparable visible light mission for hotter F, G, K stars 
for the reasons mentioned above. 

It is with these scientific goals in mind — an all sky 
survey to identify transiting planets with an emphasis 
on cool K and M stars and follow-up spectroscopic ob- 
servations using stable space telescopes — that we have 
initiated a series of laboratory experiments designed to 
probe the ultimate limits of precision photometry using 
near-IR detectors. We report on the photometric proper- 
ties of Teledyne's widely-used, infrared H2RG detectoi^ 
which of fers near unity qua ntum efficiency in the near 
infrared (Bel etic et al.||2QQ8| ). In this paper we describe 
laboratory experiments (§2) to demonstrate the techni- 
cal readiness of space-borne, transit experiments using 
existing near-IR detectors. In a future paper we will de- 
scribe the results of simulated differential spectroscopy 
planned for the JWST and proposed for other missions. 

1.1. Photometry with Infrared Focal Planes 

There is broad interest in the performance characteris- 
tics of the next generation of IR detectors. H2RGs will be 
flown on mul tiple James Webb Space Telescope (JWST) 
instruments (|Clampin|20T0l|D eming et al. 2009) and may 
also form the basis lor dedicated exoplanet observatories 
such as the candidate Explorer mission FINESSE and 
the proposed Europ ean Exoplanet Characterization Ob 



servat ory (ECHO) flTinetti et al. ||2011 



Tessenyi et al. 



12012]). Although various investigators have measured 
relevant devi ce-level properties such as the subpixel non- 
uniformity (|Barron et al. 2007 Hardy et al.|[2008), the 



count-rate dependent non -linearity or reciprocity failure 
(jBies iadzinski et al.|2011 ), and image persistence ( |Smith| 
et al.|2008[ ), there has not been a systematic study of the 
limits to high precision time-series measurements. 

Our goal is to assess the precision possible with the 
current generation of NIR devices, and whether it is pos- 
sible to achieve IR photometry that is significantly better 
than what is being achieved from space today. The cur- 
rent state-of-the-art is defined by the Spitzer and Hub- 
ble Space Telescopes. Observations to determine atmo- 
spheric composition by detecting variations in the appar- 
ent size of the planet as a function of wavelength using 
filter photometry or low dispersion spectra have had no- 
table successes, but also some controversy. Claims have 
been made about detection of wat er in Spitzer transit 
and secondary eclipse photometry ( Barmanl^OOT , 2008 




Figure 2. The substrate-removed H2RG ( |Beletic et al.|2008| > de- 
tector used in these experiments shows high quantum efficiency, 
>85%, over a broad visible and near-IR band out to 1.7/xm at an 
operating temperature of 145K. The horizontal red line show a typ- 
ical mission requirement which should be readily achievable with 
existing technology. Vertical dashed lines show a nominal system 
passband designed to limit variability due to stellar Ha emission at 
the short wavelength end (A > 0.66 /xm) and temperature induced 
detector QE variations at the long wavelength end (A < 1.65 /im). 



Tine tti et al. 2007 ) as well as counter-claims ( |Desert et al. 
2009p — sometimes using the same data set. Similarly, 
the detection of molecular features due to methane, wa- 
ter and carbon dioxide in NICMOS spectrophotometry 



is controversial (Swain et al. 2009 Gibson et al. 12011 ) 



Observations of the hot INeptune GJ 436 b have also led 
to controversial claims regarding t he presence or absence 
of methane in the atmosphere ([S tevenso n et al.| [2010 ; 
Beaulieu et al. 2011| Knutson et al. ]201l[ ). To a great 
extent it is the treatment of small, systematic detector 
effects that lie at the heart of these debates. 

The results to date suggest that the instrumental pre- 
cision of Spitzer IRAC photometry and HST NICMOS 
spectroscopy is ~ 100 ppm, although recent improve- 
ments in Spitzer's pointing stability have im proved that 
IRAC's p hotometric precision to <50 ppm (Demory et 
al. 2012). A combination of factors, including detector 
non-linearities, poor sampling of the Point Spread func- 
tion (PSF), short term wobbles and long term drifts in 
telescope pointing change the measured signal from the 
star by several percent over the course of an hour and 
make the extraction of millimag features from IRAC light 
curves very challenging. Similarly, Hubble's low Earth or- 



(Swain et al.||2008[ |2009 


Gibson et al.||2011 


Berta et al. 


2012|). 



2 H2RG or HAWAII-2RG: HgCdTe Astronomical Wide Area In- 
frared Imager 



1.2. Lightcurve Measurements 

When measuring transit lightcurves, the timescale 
of relevance is of order the transit duration, t&. An 
Earth- Sun transit produces a characteristic signal of 

depth 84^-^^^ ppm with a maximum duration 

of 13 hr. More generally, the maximum duration 
with a star-planet separation of a (in AU) is td — 
13(a/A/7) 1 / 2 (M*/M ) 1 / 2 hr assuming R* oc M*, which 
is a good approximation for the lower main sequence. 
For a solar type star transits by planets with separations 
of a ~ 0.01, 0.1 and 1 AU have durations of ~ 1, 4 and 
13 hr, respectively. Three distinct sources of noise es- 
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Figure 3. The dewar includes an LED illuminating an integrating sphere, a lenslet array, and the H2RG detector. There is no external 
source of illumination into the dewar. 



tablish the detection threshold of a transit observation 
over the measurement time t: 1) the noise from pho- 
ton arrival statistics with a standard deviation gn] 2) 
noise from intrinsic stellar variability, as] and 3) the in- 
strument measurement noise, 07. The total photometric 
noise <jt is given by 



-erf 



^2\l/2 



(1) 



In Kepler parlance gt is called th e combined differenti al 
photometric precision, or CDPP (Jenkins et al. 201 1M. 

A few general remarks are important to consider in de- 
signing a laboratory experiment to test the limits of de- 
tector performance in this application. First, the photo- 
metric noise gt is inherently a differential measurement, 
because only brightness differences in and out of transit 
are important. Thus, neither absolute photometric cal- 
ibration nor instrumental or stellar variability changes 
with timescales much longer than a few times td are rel- 
evant to the measurement. Second, over the long term, 
i.e., several transits, the stellar noise, as, should average 



down to 10 ppm f or quiescent stars (Basri et al. 2010 
Ciardi et al. 2011[ ). The Sun and some ol its G-dwar: 



cousins are relatively quiet in this respect (Basri et al. 



2010 



Jenkins et ah]p011a| ) with 10 3 ppm variation on 



timescales ol a rotation period or longer but with rela- 
tively little variability on the 10-13 hr timescale of an 
Earth transit. While some stars are as quiet as the Sun, 
10 ppm, on timescales of a few hours, Kepler has shown 
that at visible wavelengths the average solar-type star is 



somewhat noisier than the Sun, 10^30 ppm (|Jenkins et 
al.|2011b ). The vast majority of field M dwarfs observed 
by Kepler have short-term variability (0.5 to 12 hour) 
less than 1000 ppm and man y with <100 ppm (Ciardi et 



al.|2011 Bryden et al.|20lT ). However, it is important to 
note that variability due to starspots will be a factor of 
2-3 smaller in the near-IR than in the visible due to re- 
duced cojrtra£^nLti^^^y^ e ig n "J eans potion of the spec- 
trum ( Frasca et al.|2009 ). Finally, in a properly designed 
system, 07 — cr^, so that little is to be gained by reduc- 
ing any one of the noise components to be significantly 
smaller than the total. The goal of our experiments is to 
identify methods of reducing 07 so that it is comparable 
to the other irreducible terms (10-20 ppm). Quantify- 
ing the magnitude of these effects and demonstrating an 
ability to reject them as part of a total error budget is 
critical to designing a particular experiment. Major con- 
tributions to gj in transit experiments are pixel-scale de- 
tector inhomogeneities ("flat fielding errors") convolved 
with pointing changes, responsivity drifts with tempera- 
ture and other environmental effects, and detector non- 
linearities such as latent image generation. 

In §2, we describe our experiments and data acquisi- 
tion. In §3, we discuss our methods and describe the 
custom data analysis pipeline that was written to ana- 
lyze the large quantities of data that were produced by 
the experiments (up to 200 Gbytes/day). Results are 
presented §4. 

2. DESCRIPTION OF THE EXPERIMENT 
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The experimental setup described in this section was 
used to investigate three major sources of degradation 
in high precision differential photometry: 1) long term 
(~24 hr) drifts, including the response of the detector 
and warm electronics, variation in brightness of the illu- 
minating LED and environmental factors such as dewar 
temperature; 2) effects of simulated pointing drifts and 
offsets; and 3) the effects of changes in detector temper- 
ature. The magnitude of these effects and the develop- 
ment of mitigation techniques are of direct relevance to 
future space missions where the expense of controlling 
parameters such as detector temperature and spacecraft 
pointing must be traded against the exacting demands 
of transit measurements. 

A 1.7 fim cutoff HgCdTe hybrid CMOS detector with 
2044x2044 imaging pixels on 18 /am pitch, Teledyne 
H2RG-103, was mounted on a temperature-controlled 
Molybdenum plate inside a sealed chamber refrigerated 
to 77K by a liquid Nitrogen dewar with >20 hr hold time. 
The wavelength dependence of the detector responsivity 
is shown in Figure [2] The dewar was filled ~1 hour 
prior to each experiment to allow for thermal stabiliza- 
tion and thereafter was automatically refilled every 20 
hours from a 160 liter storage dewar so that the detector 
never warmed up. 

The dark chamber was evacuated through a 2.5 mm 
diameter serpentine black-painted tube to block entry of 
thermal emission from warmer parts of the dewar. Light 
leaks and thermal background were undetectable even in 
4 hour exposures against a 0.001 e~s _1 pixel -1 dark cur- 
rent floor obtained by cooling the detector to 110K. At 
140K, this detector exhibits median dark current <0.02 
e - s -1 pixel -1 which is negligible compared to all signals 
in this experiment. The ~ 20 e - /pixel read noise (for 
correlated double sampling) from the 150 pixels in the 
signal and background apertures combined contributes 
less than 1/3 of the total variance in the signal com- 
pared with the signal per pseudo star of 300,000 e - per 
1.2 sec integration. 

The layout of the dark chamber is shown in Figure [3J 
Two LEDs within the dark chamber (wavelengths 0.890 
and 1.30 jam) are driven by precision current sources ex- 
ternal to the dewar. They shine concurrently into a 2" 
integrating sphere with a finely textured gold interior and 
a precision 1 mm exit aperture. The whole illumination 
system is housed behind a baffle so that only light from 
the exit port reaches the detector 139 mm away. A sec- 
ond baffle half way along the light path blocks rays scat- 
tering off the sidewalls and shades the lenslet supports. 
A 40 mm by 40 mm fused silica lenslet array made by 
SUSS Microptics SA (Switzerland), with 300 /im lenslet 
pitch and 4.4±0.4 mm focal length was supported 5 mm 
from the detector surface to create slightly defocussed 
image of the 1 mm aperture per lenslet. The precisely 
rectilinear grid of spots was found to have a 17.29 pixel 
pitch, as expected when projection effects are taken into 
account. 

The Full Width Half Maximum (FWHM) of each spot 
was typically 2.3-2.5 pixels with low amplitude diffrac- 
tion spikes extending in the vertical and horizontal direc- 
tions (Figures [4] and [5|, as expected for a lenslet array 
with this pitch. The maximum flux per spot was ~60 
ke - s -1 corresponding to the expected count rate from a 
9th magnitude star observed using a camera with a 12.5 




Figure 4. Four of the ~10 4 spots with the scale zoomed to show 
the individual spots plus diffraction spikes. 
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Figure 5. The spot images generated by the lenslet array are 
uniform and slightly elongated. The distribution of spot sizes (Full 
Width at Half Maximum, FWHM) is narrow and centered around 
2.0 pixels in the X-direction and 2.5 pixels in the Y-direction. 

cm aperture in a very broad photometric band extending 
from 0.66 to 1.65 /im. 

A simple but precise way to move the lenslet array 
was utilized to emulate small pointing motions. Smooth, 
reproducible motion with fine resolution was produced 
by changing the gravity vector acting on the lenslet sup- 
port formed by two thin, cantilevered Aluminum sheets, 
which flex predominantly in one direction. The cylin- 
drical vacuum vessel lies horizontally in a wooden cradle 
with felt padding, allowing the orientation of the flexures 
to be changed smoothly with respect to gravity by man- 
ually rotating the vacuum vessel about the optical axis. 
A metal measuring tape attached to the vacuum vessel 
with a fixed pointer mounted on the cradle, serves as 
an angle indicator. The lenslet position is self encoding: 
while the position of any given spot can be can vary due 
to sampling effects by as much as 180 nm depending on 
registration with the pixel grid, this effect averages away 
rapidly and averages of a thousand spots exhibit 2 nm 
rms temporal stability, and comparable coherence under 
translation when comparing one half of the spots to the 
other. Figure [6] shows that the lenslet position varies as 
the sine of the indicated angle to remarkable accuracy 
for a hand driven actuator. Motion up to ±8.2 pixels is 
produced at an angle 5.2° from the row direction. The 
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0.005 pixel (90 nm) rms scatter about a sine fit is dom- 
inated by the 0.5 degree accuracy with which the dewar 
rotation angle can be read. This resolution is more than 
adequate for pointing control, while the 2 nm temporal 
stability far exceeds requirements. 

Vibration might be a concern with a flexible lenslet 
support. The resonant frequency of a mass suspended by 
a spring is uo = y/g/d, where g is acceleration due to grav- 
ity and d is displacement due to gravity. Measured peak 
displacement is 148 /im, so the resonant frequency is ex- 
pected to be 41 Hz. Indeed, the only resonance detected, 
when applying an impulse to the dewar while reading a 
small window at ~1200 Hz, was at 41 Hz. No power 
was detectable at this or any other frequency when the 
system was not excited. Nonetheless thick foam padding 
was placed between the dewar cradle and bench to in- 
crease the attenuation seventy fold at 41 Hz assuming 
moderately damped simple harmonic motion. 

The detector was read out using a 32 channel Astro- 
nomical Research Cameras controller with a 11.44 /is 
clock time, including 3 /is integration time and 1.28 /is 
settling time. Correlated double sampling (CDS) was ex- 
ecuted using conventional readout waveforms. As these 
waveforms were not optimized for such short exposures 
there was an exposure duty cycle penalty which will 
be reduced in future experiments. When idle, identical 
clocking waveforms were executed to maintain a stable 
self-heating pattern and to avoid abnormal charge ac- 
cumulation. One full frame scan was assigned to line 
by line reset, rastering through all pixels with identical 
timing to the readout so that the self heating footprint 
during reset matched the read scans. This was immedi- 
ately followed by a read scan to measure the post-reset 
offset, then an identical read scan to measure integrated 
signal plus offset. The difference of the two samples for 
each pixel is the signal. The exposure time is the delay 
plus one frame scan time, while the exposure periodicity 
is the delay plus three frame scan times. The exposure 
duty cycle is their ratio, 33%. To improve the duty cy- 
cle, while avoiding some cosmetic defects, the frame scan 
time was reduced from 1652 msec for a full frame to 1221 
msec, by skipping 400 lines and then reading a band 1513 
lines high. 

A typical observing sequence consisted of filling the 
dewar, setting the LED intensity, and starting a steady 
cadence of one 1.2 sec exposures every 3.6 sec, repeated 
continuously for periods as long as the 20 hr LN2 autofill 
interval. This experiment was repeated with the same 
timing for a range of illumination intensities. In one 
test, the intensity was deliberately modulated to mea- 
sure the effectiveness of normalization by the ensemble. 
In other cases, the dewar was rotated to move the im- 
ages, simulating pointing changes. Other disturbances 
such as temperature changes were introduced to assess 
their importance and our ability to correct for them. 

3. THE ANALYSIS: METHOD AND DATA PIPELINE 

Raw image data for each experiment were run through 
a data analysis pipeline that was written specifically for 
the experimental layout described in §2. Individual steps 
of the pipeline are summarized below. 

• Define a region of interest on the array selecting an 
area of high cosmetic quality for both the detector 



and lenslet arrays, roughly 30 by 30 spots in size. 
For analysis we typically used 1,000 spots avoid- 
ing bad pixels. Channel-to-channel variations in 
the warm electronics readouts can introduce small 
measurement perturbations which can be reduced 
by rigorous design of the warm electronics and 
standard destriping techniques. Since for these ex- 
periments we wanted to concentrate on the fun- 
damental detector performance issues, we avoided 
pixels spanning channel boundaries. 

• Construct a master dark frame from 200 continu- 
ously acquired dark frames and subtract from each 
image. No flat field correction was applied to the 
data since the analysis steps described below pro- 
vided the equivalent of a highly localized flat field 
for each spot. 

• Determine centroids for each star for use as the 
center for aperture photometry. The IDL routine 
FINE0 locates sources by fitting a 2-D Gaussian, 
yielding the centroids and full width at half maxi- 
mums (FWHMs) in orthogonal directions for each 
star. Depending on the experiment, we developed 
two different ways of placing apertures. In experi- 
ments where the stars are stationary with respect 
to the focal plane, i.e. with no intentional pointing 
disturbances, we used stationary apertures which 
remained at a fixed position throughout the pho- 
tometric sequence. To account for any drifts that 
may occur, as well as any estimation errors that 
may arise due to quantization effects of a spot ly- 
ing on a pixel boundary, the aperture centers are 
computed as follows: the first few (~ 5) frames of 
the sequence are averaged together and then FIND 
is used to locate the spots in this averaged frame. 
The same is done with frames from the middle 
and end of a sequence. This produces centroids 
for each spot which are then averaged together to 
get the final value of the aperture center. Typical 
image motions from the beginning to the end of 
the photometric sequence are ~ 0.03 pixels and ~ 
0.004 pixels for X and Y centroids, respectively. In 
experiments with deliberate pointing disturbances, 
the apertures are allowed to track the star and the 
moving centroids were recomputed every frame. 

• Calculate fluxes, using the IDL routine 
APER^Jwith a 3 pixel radius aperture which cap- 
tures ~ 95% of the point spread function's encircled 
energy. The signal contained within a background 
annulus of inner and outer radii of 5 and 8 pixels 
was subtracted to remove background counts. 

We first use aperture photometry of 10 spots to demon- 
strate that the detector performance is limited by photon 
shot-noise. Measurements at 5 different brightness levels 
were fitted to a noise model (Eqn 2) relating the mea- 
sured noise to the signal intensity and read noise (RN) 
in N P i X = 150 pixels (signal aperture plus background 
annulus): 

Noise(e-) = Signal(e~) + (RN) 2 N~ (2) 

3 IDL Astronomy User's Library 
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Figure 6. As the dewar is rotated around its horizontal axis, the average position of the spots shifts by ±8 pixels along the rows (X- 
direction). A simple sinusoidal model provides a good fit to the data with residuals at the level of a few hundredths of a pixel. Note that 
the shifts along columns (Y-direction) are plotted on a 10-times exaggerated scale. 

af over a time lag r is the second moment of L + — L~ , 
o?(t) = \((L+-L-) 2 ) Vr<|, (3) 

where L + and L~ are integrals (or equivalently, the in- 
tegrals are replaced by sums for periodically sampled, 
discrete light curves) over adjoining time intervals of du- 
ration r, 




Figure 7. Measured noise in 10 photometric apertures is shown as 
a function of signal brightness along with a a simple model combin- 
ing read noise and photon shot noise (Eqn. 2). The measurements 
are photon shot noise limited. 

where the Noise and Signal in Data Numbers (DN) were 
converted to e~ using a conversion gain, G. Dark cur- 
rent is negligible at these signal levels. The solid line in 
Figure [7] shows the noise model with a conversion gain of 
G = 4.2 ± 0.2 e~/DN and RN = 15 ± 2 e~ . The conver- 
sion gain includes a correct ion for interpixel capacitance 
of 0.81 ( [Brown et al~|[2Q06 ). 



- / l(t)dt 

T Jt-T 

u 



l(t)dt 



and T is the total duration of the lightcurve. This statis- 
tic preferentially selects noise characteristics on a lag 
timescale, ~ r, while rejecting noise on timescales both 
longer and shorter than r. For random noise, the vari- 
ance assumes its usual 1/r dependence; and in general, 
f~ n noise has af ex r n_1 . It is illustrative to write the 
Fourier domain filter-form of the Allan deviation, as 



3.1. Light Curve Statistics 



PCX) 

2 / ^(/)sinc 2 (^r/)sin 2 (^r/)d/ 
Jo 



We use the Allan variance (Howe et al. 



1981D as 

the main measurement statistic tor assessing photomet- 
ric precision. For a light curve l(t), the Allan variance 



where Si(f) is the power spectrum of l(t). This shows 
that Si(f) is filtered by a bandpass filter with a response 
H A (f), 
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\HaU)? 



) shv^7rr/) 



matrix with elements given by the partial derivatives, 



peaking at / c± 1.2/ttt. 

For a given experiment, we calculate an Allan devi- 
ation for each star within the detector region of inter- 
est. The averaged or global Allan deviation is used as 
the statistic for the determining how well we did over 
the ensemble of stars in the experiment and is given as, 
&l(t) — \J (of (t)), where the average () is over the en- 
semble, after rejection of outliers that show exceptionally 
poor Allan deviations. About ~ 3 — 5 % of all stars are 
rejected as outliers with anomalous Allan deviations. In 
this paper we do not discuss the nature of the the outlier 
stars other than to say that their spatial distribution is 
fairly random and that at least some fraction of these 
stars show poor photometry due to unresponsive pixels 
located within the aperture. 

When no corrections are applied to light curves, cr L (r) 
quickly departs from the r -1 / 2 curve expected for Gaus- 
sian noise. The primary culprit at this noise level are 
the ~ 10 3 ppm fluctuations of the lamp on timescales of 
minutes (Figure |8k. Because the lamp intensity is not ex- 
plicitly stabilizeaby feedback, we rely on common mode 
rejection of lamp intensity variations as described below 
to reduce this noise source. 

3.2. Measurement Model and Parameter Decorrelation 

Light curves rapidly depart from the ideal case due to 
a number of sources of 1/ f noise. Laboratory lamps vary 
more rapidly than stars, and lamp fluctuations (10 3 ppm) 
are typically the largest source of non-ideal behavior in 
our experiments. However, because the lamp intensity is 
in common to the spot ensemble and because we achieve 
excellent common-mode rejection, lamp variations can 
be removed to the required accuracy. The stability of 
the light curves is also affected by environmental factors 
such as image wander, changes in the optical point spread 
function, detector temperature, etc. These changes can 
be measured and captured in the "instrument state" and 
then be removed from the data. In general, the light 
curves y (vector-valued for m spots) are related to the 
combined source and instrument state, x, as 



y = F(x). 



(4) 



There are various ways in which to define the vector- 
valued model function F(x). The most sophisticated ap- 
proach is via a full description of the image formation 
process if one is in possession of the detailed PSF and 
the detector response function. Such an approach may 
be necessary to obtain good photometric performance in 
cases when the PSF is poorly sampled and/or if the per- 
turbations are relatively large. In our case, with a well 
sampled PSF, and for relatively small perturbations to 
the instrument state, a simpler parameterized form of 
F(x) suffices. We find that in most instances we can 
linearize (3) and rewrite it as 

y = F(x i )+K(x-x i ), (5) 
where x, is some initial state, and K is a m x n Jacobian 



Ki. 



dxi 



(6) 



relating the photometry for spot i to the instrumental 
parameter j. One of our goals is to quantify the "sen- 
sitivities" to various parameters; for any parameter Xj, 
this is simply a row average (i.e. equivalent to an average 
over the spot ensemble) of the Jacobian, (Kij)i. 

However, in some cases the truncation to first order 
(4) is insufficient; this is the case if there is relatively 
large image motion (> 0.1 pixel) or if lamp fluctuations 
are large and pixel response non-linearities must be con- 
sidered. If so, F can be expanded to second order in x. 
This higher order model is recommended, for example, 
for light curves extracted fro m poorly sampled Sp itzer 
IRAC Bands 1 and 2 images (Knutson et al. 2011), re- 
quiring estimation of the Hessian, 



Hi 



dxjdxk 



(7) 



In our experiments the state is well determined by the 
following set of parameters: 

1. The image centroid or pointing drift, 5x and 5y, 
measured in pixels, and usually a small fraction of 
the pixel pitch. In general, we find that for image 
motions larger than a tenth of a pixel, 5x,5y > 0.1 
pixels, second order terms in Sx and Sy are required 
to model the photometry. With image motions of > 
0.3 pixels the second order model no longer suffices. 

2. Changes in the optical point spread function. Fo- 
cus drifts can lead to significant changes in the PSF 
and the photometry. We capture focus with rel- 
atively simple parameterization, i.e., we approxi- 
mate all spots as Gaussians and track the change 
in the FWHM (5w x , Sw y ) caused by temporal fo- 
cus drifts. However, we have no way to inject con- 
trolled amounts of focus error. 

3. Changes in the lamp intensity (5L) are estimated 
by an average of the lightcurves of an ensemble of 
stars randomly distributed over the detector ROI. 
In typical experiments, about 20 stars are need for 
a reasonable estimate of 5L. 

4. Any changes in the detector temperature (STfpa)- 

Often changes in the state parameters are correlated, 
being driven by the same underlying mechanism, i.e., the 
daily thermal cycles between dewar fills. On orbit this 
could be the solar thermal forcing of a spacecraft. Given 
these correlations it is convenient to transform the state 
x to x 7 , with the latter having a minimum set of mi- 
correlated or orthogonal parameters by using a principal 
components analysis (PC A). PC A is an orthogonal linear 
transformation that maps a set of observations of possi- 
bly correlated vectors into a set of values of orthogonal 
vectors, called the principal components. By definition 
the first principal component has the highest variance, 
i.e. it accounts for as much of the variability in the ob- 
served data as possible. Subsequent principal compo- 
nents are ordered by variance but under the constraint 
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4000 



Figure 8. The light curve of a typical spot for the first hour after turning on the LED shows a rapid drop in LED brightness of ~1.2%, 
as expected for diode junction self heating by ~1K, followed by a long term drift of — 0.2%/hr. Short term fluctuations at the level of a 
few parts per thousand are also evident. 



of orthogonality to the preceding components. Our anal- 
ysis proceeded according to following steps: 

• Let the light curve for each spot consist of p mea- 
surements; each measurement yields an estimate of 
the state parameters. Convert the p x 1 vectors of 
state parameters Xj to centered, rescaled versions 



via the transformation 



(xj - Xj)/cr Xj , where 



Xj and a Xj are the mean and standard deviation. 

• Construct the p x n state matrix X for each spot 
using a set of n state parameter vectors. Rows of 
X correspond to the observations, and columns to 
the parameters. 

• Use the IDL routine PCOMP(X)F]to perform prin- 
cipal components analysis on a and return the 
principal component coefficients. 

• We generate the new state, X, using a truncated 
set of m principal component vectors (m < n, and 
m is typically 3-4). The normalized light curve Z is 
described by Z = XK + e, where e is an error vec- 
tor. Minimizing the cost, e T e, with respect to K, 
yields a least-squares solution for the sensitivities 
K = (X T X)- 1 X T Z. Repeated tests show the er- 
ror e to be nearly zero mean Gaussian, making the 
least squares estimate of K equivalent to a maxi- 
mum likelihood estimate, and conferring upon it a 
direct probabilistic interpretation. 



• Compute the model light curves, XK. The error 
residuals, e = Z — XK are used to compute the the 
Allan deviations. Allan deviation curves labeled as 
"decorrelated" have been corrected for systematics 
using the above scheme. 

4. RESULTS 

In the simplest experiments, we held the experimen- 
tal setup constant for approximately 24 hr and either 
kept the lamp brightness constant or changed the lamp 
brightness in discrete steps. Figure [9^i shows the inten- 
sity of a single spot before correction for variations in 
lamp brightness. The Allan deviation as a function of 
lag for that single spot (Figure |9fc>), inclusive of lamp 
variations, hits a floor at the 0.1 to (or 1000 ppm) level. 
The decorrelation using PCA (with or without clipping of 
a few bad pixels) greatly reduces the effect of lamp vari- 
ations and other noise sources, achieving a noise floor of 
~50 ppm for a single spot. Data sets from different ex- 
perimental runs shown in Figures To][l3 show single spot 
noise floors between ~30-50 ppm. 

Figure [T0| shows the result of first averaging the light 
curves oirnultiple spots together and then computing 
the Allan deviation of the averaged light curve. At short 
lag times, the initial noise level decreases as expected 
for photon-noise limited measurements: the noise is in- 
versely proportional to the square root of the number of 
spots, i.e. proportional to the square root of total number 
of photons. At longer lag times, the limiting noise floor 
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also drops as the square root of number of spots, i.e. as 
the square root of number of independent, uncorrelated 
measurements. This result implies that averaging over 
multiple spots can reduce other temporally-correlated 
noise sources down to the ultimate noise floor seen in 
these experiments of 10-20 ppm. Figure [TT] shows that 
even large, commanded excursions in lamp brightness 
can be corrected for via ensemble averaging and PCA 
decorrelation. Overall, these data show that by using 
decorrelation techniques and suitable averaging of inde- 
pendent samples, the H2RG detector can achieve very 
high precision differential photometry with noise values 
as low as £ 20 ppm. 

4.1. Pointing Rejection 

In general, the measured count rate for a source repre- 
sents a convolution of the detector spatial response with 
the pattern of optical illumination over a predefined aper- 
ture. The illumination pattern is the broadband PSF in 
the case of a survey experiment, or a narrow spectral ex- 
traction region in a spectrograph. Because of the convo- 
lution the coupling of pointing changes and/or PSF vari- 
ations with inter-pixel or intra-pixel responsivity varia- 
tions within the detector changes the measured signal. 
While a transit survey experiment can correct for gross 
drifts in detector /electronics response using stellar en- 
semble averages, changes in the measured brightness of 
individual stars due to small pointing shifts must be cor- 
rected using decorrelation techniques. These techn iques 
have bee n pioneered using Spitzer and HST data ( |Knut-| 
son et al.||2007||Swain et al.|2QQ8| >. 

The response within a pixel can have substantial struc- 
ture due to declining detection efficiency at its edges, 
diffusion of charge into neighboring pixels and electrical 
crosstalk between pixels. Sensitivity variations across a 
pixel reduce photometric accuracy in images even where 
the PSF is fairly well sampled. The intrapixel variations 
in the PACE type detect ors flown on N ICMOS have been 
found to be quite large ( |Finger et al.||2000| ) and clearly 
affect the lightcurve fidelity obtained during NICMOS 
transit photometry, despite the large (5-pixel FWHM) 
defocus in the NIC3 camera. 

These intra- and inter-pixel responsivity variations 
make pointing jitter and wander amongst the largest 
terms in the photometry error budget in a transit ex- 
periments using older detectors. Improvements to quan- 
tum efficiency of HgCdTe detectors in recent years have 
resulted in most photons being detected if they are not 
reflected. Even though lateral charge diffusion results 
in a variation in response of a given pixel depending on 
the position of the stellar image the integrated response 
remains almost constant. We conducted an experiment 
to determine the degree to which pointing affects photo- 
metric precision, and the degree to which decorrelation 
analysis can be used to remove these effects from the 
lightcurves. Controlled pointing offsets were introduced 
by rotating the dewar and causing a 1-D shift in spot 
positions (§2, Figure [6]). A series of dewar rotations pro- 
duced up to 0.25 pixel steps in the X locations of the 
spots and much smaller shifts in the Y-locations. After 
determining the new spot centroids, we applied the PCA 
decorrelation technique and were able to reduce the noise 



floor to well under 30 ppm as demonstrated in Figure [12 
We conclude that in a configuration with ~ 2 pixel ra- 



dius spots, photometry of <50 ppm can be achieved using 
PCA correction even in the presence of pointing offsets 
of up to 0.25 pixel. 

4.2. Temperature Dependence 

Thermal fluctuations can produce changes in detector 
response. Fortunately, the temperature dependence of 
detector responsivity is relatively weak in modern H2RG 
arrays compared with variation in the quantum efficiency 
in the NICMOS 3 devices used on HST. We conducted a 
simple experiment to determine the degree to which the 
detector temperature must be controlled to avoid impact- 
ing precision photometric requirements. In our set-up a 
temperature controller maintains the focal plane at 140 
K with a precision of ±5 mK. An 8.5 hr long photometric 
image sequence was acquired during which the temper- 
ature of the focal plane was modulated by changing the 
controller set-point temperature by up to 2 K. The tem- 
perature steps are shown in Figure [13] along with the 
change in the average signal. The figure also shows the 
Allan deviation for the uncorrected and corrected (decor- 
related) data. The decorrelation technique can correct 
for relatively large swings in detector temperature with- 
out a significant impact on the final photometric preci- 
sion. 

The average temperature coefficient of the intensity 
variation found by averaging the lightcurves of many 
stars is ~ 250 ppm/K. The combination of temperature 
control at the 0.1 K level and decorrelation analysis will 
reduce the effect of temperature fluctuations to accept- 
able levels. A flight instrument would use a variety of 
techniques to maintain detector temperature including 
heaters at the focal plane level (0.05 K) as well as us- 
ing thin metal traces built into the H2RG multiplexor 
that can monitor and control the package temperature 
directly at the few milli-K level. Finally, detector refer- 
ence pixels can be used to suppress offsets due to residual 
bias voltage and temperature drifts. 

5. DISCUSSION 

5.1. Applicability to Other Near-IR Focal Planes 

How applicable are these results to Teledyne HxRGs 
with different cutoff wavelengths? H2RG-103, used in 
these tests, is from one of the first successful batches of 
1.7 /im cutoff HgCdTe detectors made for the SNAP pro- 
gram in 2006. Substrate removal had been mastered, de- 
livering what is believed to be near 100% internal quan- 
tum efficiency: almost all photons passing through the 
anti-reflection coating produce signal electrons. A corol- 
lary of this high internal QE, shared by devices with 
other cutoff wavelengths, is that intrapixel response vari- 
ation is low since there is no charge loss mechanism to 
modulate sensitivity. 

Barron et al. (2007) performed intra-pixel response 
measurements on the near-identical 1.7 /am cutoff H2RG- 
102 from the same batch, and found that, while lateral 
charge diffusion redistributes charge between neighbor- 
ing pixels, the integrated flux is independent of the posi- 
tion of sub pixel size spots at the <2% rms level. Hardy 
et al. (2008) reached the same conclusion after perform- 
ing; similar tests over a wider area on JWST-like HIRG^ 

4 H1RG or HAWAII-1RG: HgCdTe Astronomical Wide Area In- 
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Figure 9. a, left) A single spot light curve (h/I — 1) over an 8 hr interval for nominally constant lamp intensity shows variation at the 
level of a few tenths of a percent (black curve). The PC A model is shown in magenta, b, right) The Allan deviation of the light curve 
(black trace), inclusive of lamp variations, hits a floor at the 0.1% (or 1000 ppm) level. After PC A decorrelation the noise level obeys a 
r -0.5 relationship until reaching a floor of roughly 50 ppm after 4.2 hours for single spot before (blue curve) and after bad pixel rejection 
(green curve) 
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Figure 10. Allan deviations of decorrelated light curves, based on data presented in Figure [qJl made by co-adding multiple spots (from 
top to bottom: N = 1, 9, 25, and 100 spots). The noise level decreases as \/N at short lags where the statistics are determined by photon 
noise. Even at longer lags, where correlated noise dominates single spot Allan deviations, the improvement is ~ y/N, suggesting that the 
systematics are mostly uncorrelated between spots. The dashed line shows that when sufficient flux is included in the measurement, i.e. 
the 100 spot line, the noise maintains the photon noise limited r -0,5 relationship until reaching a floor of roughly ~ 10 — 20 ppm. 

( Beletic et al.|2008 ) with 5 /im cutoff and found that pro- sent a worst-case scenario, and that the results in this 
cess dependent features, such as surface particulates and paper are probably extensible to longer cutoff devices. It 
sub pixel width scratches dominated the response. is often the case that greater differences in performance 



The higher bandgap of the 1.7 /im detector presents 
special manufacturing challenges that result in higher 
read noise and persistence. While operating tempera- 
ture can be higher for a given dark current, the dark 
current floor is similar to the longer cutoff devices. For 
these reasons we believe that the 1.7 /im devices repre- 

frared Imager 



are observed due to process variations than due to the 
cutoff itself, but the better devices at any cutoff wave- 
lengths are expected to reach the levels of performance 
reported here. 

Thus, we maintain that our results showing 1.7 /im 
cutoff H2RG devices are capable of < 50 ppm precision 
photometry are directly applicable to longer wavelength 
devices such as JWST's 5 jam arrays. Our existing ex- 
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Figure 11. left) The light curve of a single spot following commanded changes in lamp intensity, right) The averaging characteristics of 
unnormalized data are in black with 'x' markers. After normalization by the ensemble average, the noise floor drops to <100 ppm (blue 
curve, diamonds). After full PCA decorrelation, the noise level follows a r~ 0-5 relationship until reaching a floor of roughly <50 ppm for 
a single spot (green curve, triangles). The common mode rejection for lamp variations is better than 100:1. 




5,000 10,000 15,000 20,000 25,000 30,000 TjlTlB (S6C) 

Figure 12. left) Light curve of a single spot in response to commanded image motion. The large steps (right hand axis in pixel units) 
represent centroid motions of up to 0.25 pixel. The thick black curve shows the resultant variations in measured signal intensity, up to 
almost 1%. The lighter colored (blue) curve shows the PCA model fitted to the light curve data, right) The curves marked with squares 
denote data with no motion while the curves with triangles denote data with position offsets. The two lower curves were corrected using 
PCA analysis. After full PCA decorrelation the noise level follows a r~ 0-5 relationship until reaching a floor of roughly 30 ppm for a single 
spot. 




Figure 13. left) A single spot light curve responding to commanded steps in the detector temperature in black, right) An Allan deviation 
plot showing unnormalized data (black curve) at few percent level. After normalization by the ensemble average, the noise floor drops to 
<100 ppm (blue curve). After full PCA decorrelation the noise level follows a r -0,5 relationship until reaching a floor of roughly <30 ppm. 
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perimental setup could be used to probe the behavior of 
2-2.5 /im cutoff material with only slight modifications. 
However, operation out to 5 /im would require a new ex- 
periment designed with much more attention to reduced 
backgrounds and lower temperature operation. 

5.2. A Near-IR Transit Survey 

We put this detector performance into context by con- 
sidering observations with a small near-IR camera (12.5 
cm diameter) observing in a wide "Jy^" band extend- 
ing from 0.66 to 1.65 /im where we have taken typical 
values for the overall optical and detector quantum effi- 
ciency (total of 50%) and sky brightness in space (0.15 
MJy sr _1 ). In Figure [l4| the solid black curve shows 
the idealized case wherethe SNR varies simply as the 
square root of the stellar photon flux. Three additional 
curves show realistic cases taking into account the pho- 
ton noise from sky background and 3 different values for 
the noise floor (10, 50 and 1,000 ppm). For the bright 
stars considered here, detector read noise and dark cur- 
rent are negligible sources of noise. For the brightest 
stars, the SNR curves flatten out due to the uncorrected 
detector responsivity changes, i.e. the flat field "noise 
floor". At faint magnitudes, photon noise from the sky 
background reduces the SNR relative to the stellar pho- 
ton noise limit. The 1,000 ppm limit represents the best 
values that have been obtained in ground based surveys 
of very large numbers of stars (§1). Horizontal lines de- 
note the SNR levels needed to detect the transit of a 1 
R e planet orbiting G2, K5, and M3 stars at the 3<r level. 
The noise floor demonstrated in this paper, <50 ppm is 
adequate to detect at Earth transiting a K star at 10-20 
pc in a single 900 sec observation where 900 seconds is 
a typical observation time for a spacecraft in Low Earth 
Orbit. The fact that the noise in the average of indepen- 
dent measurements continues to decrease to ~10-20 ppm 
(§4) implies that weaker transit signals can be extracted 
from the light curves with suitable processing of multi- 
ple transits. Typical requirements for a transit survey 
demand a minimum of 3 transits and a total SNR=7 for 
a confirmed detection. 

The yield of a survey depends on many instrumental 
and survey strategy parameters, e.g. number of cameras, 
survey cadence and duration, mission length, etc., and 
is beyond the scope of this article. We have developed 
a Monte Carlo simulation using modest extrapolations 
of the distributions of the size and orbital location of 
planets derived from radial velocity and t ransit surveys 



(Howard et al. 2011 



Borucki et al. 2011). The simula- 



tion, which will be described elsewhere, conservatively 
estimates that for a sample of 2 x 10 6 stars each hosting 
exactly one planet and spanning spectral types from F0 
to M4 (J w < 14 mag) and M5-M8 (J w < 15 mag), a 
1.5 yr near-IR survey using 6 cameras on a spacecraft in 
a low earth orbit would detect over ~200 rocky planets 
( ^ 2 R®) (Figure 15) of which at least ~50 would be lo- 
cated in the habitable zones of their K and M star hosts. 
The time any star is monitored varies with ecliptic lat- 
itude, /3, and for the particular configuration simulated 
here we assumed 75 days per year for a star at j3 = 30°. 
The survey would also identify a few thousand gas and ice 
giant planets covering a broad range of equilibrium tem- 
peratures. Some 300 planets would orbit stars brighter 
than Jw = 9 mag and be suitable for confirmation and 



mass determination via radial velocity measurements and 
spectroscopy for the characterization of planetary atmo- 
spheres. A survey conducted in a more distant orbit from 
the Earth, e.g. L2, would offer a greatly improved ob- 
serving cadence resulting a more complete survey and 
better characterized transit curves. 

Due to the fundamental selection biases of transit sur- 
veys of modest duration, most transiting planets will be 
found in short period orbits of a few months or less, but 
only for late type stars does this preferred orbital loca- 
tion coincide with the habitable zone. A near-IR survey 
(§1, Figure fl| would find 3-7 times more rocky plan- 
ets orbiting Tv and M stars than an otherwise identical 
visible light survey, a critical difference for the identifi- 
cation of planets with temperate atmospheres (250 K < 
T e ff < 350 K) which are of great interest for spectro- 
scopic follow-up. The laboratory experiments described 
in this paper demonstrate that the performance of mod- 
ern near-IR detectors is adequate for such a survey and 
for subsequent follow-up observations. 

5.3. Future Work 

Experiments with the existing dewar configuration are 
currently underway to correlate spatial variations detec- 
tor performance (e.g. dark current, responsivity, and lin- 
earity) with photometric performance, both to identify 
whether calibrations can improve upon the photometric 
performance obtained through PCA alone and to estab- 
lish criteria for excluding "bad pixels". These criteria 
will assist in developing metrics for detector screening 
and grading for possible future missions. Future experi- 
ments in this series will include 1) an examination of the 
effect on photometric precision of varying the number 
of pixels within the stellar image which is directly rele- 
vant to JWST/NIRCam observations using weak lenses 
to defocus NIRCam images and 2) simulations of the ex- 
traction of transit spectra over a wide range of spectral 
line properties, e.g. line depth and line width. These 
results will be relevant to photometric, grism and spec- 
troscopic modes of JWST as well as to dedicated transit 
spectroscopy missions presently under consideration in 
the US and Europe. 
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Figure 14. Signal to Noise ratio (SNR) for stars as a function of broad band (0.66-1.70 /im, denoted J\y) magnitude after a 900 sec 
integration typical of an observation in Low Earth Orbit. Observations are assumed to use a camera with 12.5 cm aperture using a 1.7 
fim cutoff H2RG detector as discussed herein. The noise model includes photon noise from the star and the sky as well as a noise floor 
set by the limiting detector performance (the "Flat Field" limit, FF) in ppm. The solid black line denotes the noise limit as set only by 
stellar photon noise. At magnitudes fainter than Jw ~14 mag sky background reduces the SNR while for bright stars the residual noise 
floor dominates the SNR. The effects of stellar jitter and background confusion are neglected here but are included in the Monte Carlo 
simulations. The horizontal lines show the SNR needed to detect a 1 Rq planet transiting a G2, K5, M3 star in a single 900 sec observation 
at the 3cr level. 
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Figure 15. A Monte Carlo simulation suggests that an all sky survey would find thousands of gas giant planets and hundreds of rocky 
planets, many of which would occupy the temperate habitable zones of their late-type host stars. In this polar diagram, the radial coordinate 
represents a J magnitude and the polar angle represents a simulated sky coordinate. Gas giants from the simulation are shown as small 
filled circles while rocky plane ts are shown as red 0s. For comparison, confirmed Kepler planets and Kepler candidates are shown with 
their J (2MASS) magnitudes ( |Borucki et al. 2011 Batalh a et al. |[2012| . Kepler's gas giant planets are shown as blue circles and rocky 
planets as red squares. 
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